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Abstract 

The space subdivision in cells resulting from a process of random nucleation and growth is a 
subject of interest in many scientific fields. In this paper, we deduce the expected value and variance 
of these distributions while assuming that the space subdivision process is in accordance with the 
premises of the Kolmogorov-Johnson-Mehl-Avrami model. We have not imposed restrictions on 
the time dependency of nucleation and growth rates. We have also developed an approximate 
analytical cell-size probability density function. Finally, we have applied our approach to the 
distributions resulting from solid phase crystallization under isochronal heating conditions. 
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I. INTRODUCTION 



In this paper we consider the subdivision of a D-dimensional Euchdean space into disjoint 
regions created after a process of random nucleation and growth. Random subdivisions can 
be obtained by several different methods, amongst which Poisson-Voronoi and Johnson- 
Mehl tessellations^ have been widely studied. The Poisson-Voronoi tessellation is obtained 
by randomly picking several points, the seeds Pi, by a Poisson process. Next, the space 
is subdivided in cells, Ci, by the rule: Ci contains all points in space closer to Pi than to 
any other seed. This cellular structure is extensively applied in many diverse scientific fields 
including biology,-i^ computer science,-i^ materials science,-*^ astrophysics,->^ii^ medicine,— 
agriculture,— quantum field theory,-- and sociology.— 

The space tessellation can be fully characterized by means of the probability density 
function (PDF), f{s), which is the probability that a cell has a size between s and s + ds. 
The properties of the PDF of the Poisson-Voronoi tessellation have been extensively studied 
both theoreticallyi'^i^ and numerically.-'i^'^'^i^'22i2ii2I It is well known that the Poisson- 
Voronoi tessellation PDF is described by a gamma distribution 



where F is the gamma function, u is a parameter that is dependent on the dimension D, i.e. 
= 2, 3.584 and 5.586 for D = 1, 2 and 3 respectively, and E is the expected cell size. 



It is worth mentioning that Eq. ([T]) has been analytically derived for the one-dimensional 
case where 1^=2 is an exact result. Conversely, for the two and three-dimensional cases, 
the validity of Eq. ([1]) is supported by analytical approximations and numerical fits. 

Our main interest is the characterization of grain morphology related to crystallization. 
In general, the crystallization of most materials takes place by means of a nucleation and 
growth mechanism: nucleation starts with the formation of small atom clusters of the new 
stable phase in the metastable phase. Subsequently, clusters with sizes greater than the 
critical, or nuclei, start to grow by incorporating neighboring atoms of the metastable phase. 
During this growth, grains impinge upon each other. Finally, the structure of the new stable 
phase consists of disjoint regions or crystals separated by grain boundaries. The evolution 
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of crystallization and grain size distributions is entirely determined by the nucleation rate 
density / and the grain linear growth rate G. When nucleation takes place for a very short 
time, its rate may vanish before the onset of particle growth {site saturated nucleation) 
In this case, the crystal structure is equivalent to a Poisson-Voronoi tessellation provided 
that nucleation is Poissonian through the whole space and growth is isotropic. 

Conversely, continuous nucleation takes place when nucleation and growth occur at the 
same time. In general, there is an energy barrier for nucleation and growth to happen. 
Thus, I and G depend on temperature. For the particular case of isotropic and isothermal 
transformations, where / and G are constant, the resulting crystal structure corresponds to 
the well-known Johnson-Mehl tessellation."! For this tessellation. Axe et al.— have obtained 
an analytical solution for the one-dimensional case while Mulheran^ has developed a simple 
(but not so accurate) relation for the two- and three-dimensional cases. Alternatively, Monte- 
Carlo simulations provide a powerful tool for the calculation of tessellations and PDFs under 
a wide variety of conditions.-i^Ii^^i^iSSiSii^Si^^iMi^^ 

Under non-isothermal conditions, / and G depend on time by virtue of their tempera- 
ture dependence. Therefore, an infinite number of different tessellations/structures can be 
obtained by varying the thermal history. Unlike the Poisson-Voronoi and Johnson-Mehl 
tessellations, the analytical results related to tessellations emerging from time dependent 
nucleation and growth rates are scarce. Indeed, as far as we know, the analytical models are 
limited to time dependent nucleation rates.— In particular, Jun et al.™ have derived an 
analytical solution for the one-dimensional case. Particularly relevant to the present work 
are the results of Pineda et al.— who have obtained an accurate analytical description for 
the two- and three-dimensional cases. 

In the present work we will consider those transformations that fulfill the Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) premises. No restrictions will be imposed on nucleation 
and growth rate time dependence. We will refer to these tessellations as KJMA tessel- 
lations. KJMA theory has been widely applied to describe systems undergoing first-order 
phase transformations. For instance, DNA replication;^^ crystallization of polymers,— amor- 
phous materials^Ss^ and glasses;^ switching in ferroelectrics^^ and ferromagnets;^^ lattice- 
gas models;— and film growth on solid substrates^^i In Section II we will describe the basic 
concepts of KJMA theory and will focus our attention on those aspects that are useful to 
the development of our work. Section III is devoted to the calculation of the expected value 
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and variance of the distributions related to the KJMA tessellations. In Section IV we will 
derive a simpler approximate relation for the variance and will check its accuracy. As an 
application of the previous results, in Section V we will derive an approximate grain size 
PDF which is the superposition of gamma distributions. Finally, at the end this section 
we will verify that the grain radius PDF can be expressed as well as the superposition of 
Gaussian distributions. 

II. THE KOLMOGOROV-JOHNSON-MEHL-AVRAMI THEORY 

The KJMA theory^"^"^ describes in a very simple form the kinetics of transformations 
governed by nucleation and growth that satisfy the following assumptions: 

i nucleation must be Poissonian through the entire space; 

ii the volume of an arbitrary grain is much smaller that the volume of the system; 

iii the crystal growth rate is isotropic. 

On the basis of these premises, Kolmogorov calculated the evolution of the transformed 
fraction, X{t), through the probability, p{t), that an arbitrary point O has not crystallized, 
i.e., the probability that no nuclei able to transform O will be formed during the time interval 

M, 

Xit) = l-p{t), (3a) 
-go [ I{r)r{t,TfdT 



p{t) = exp 







(3b) 



r(t,r) = j G{z)dz, (3c) 

where gr, is a geometrical factor related to the shape of the crystal - for a D-dimensional 
sphere go = vr^/^/r(D/2 + l) - and r{t, r) is the minimum distance between O and a nucleus 
created at r, so that the nucleus would not transform O. 

Based on geometrical arguments, Avrami deduced the following relation: 

dMt,T) ^ 1-X{t) 

dtv^t^r) l-X(r)' ^> 

where dtv{t,T) and dtVex{t,T) are respectively the actual and extended average volumetric 
growth rate at time t for grains nucleated at time r. The word extended refers to the 



volume a grain would attain if nuclei grew through each other and overlapped without 
mutual interference. 

The integration of Eq. (jlj) leads to^ 

^-"X.At). (5) 

Finally, integration of Eq. ([5]) gives Avrami's well-known formula 

X{t) = l-exp[-Xe.(t)]. (6) 

The calculation of X^xit) is straightforward and obtained by simply neglecting the impinge- 
ment between nuclei 

Xex{t)=gD [ I{r)r{t,TfdT. (7) 
Jo 

The combination of Eqs. ([6]) and ([7]) gives Eq. ([3]). As it is well known, Avrami and Kol- 
mogorov deduced the same relation using different approaches. 

Note that in Eq. ([7j) it is assumed that the nucleation rate is not affected by the shrinking 
of the untransformed phase. In the calculation of Xex{t) the phantom nuclei are taken 
into account. Avrami designated as phantom nuclei those nuclei that are formed in the 
transformed fraction and therefore do not contribute to the formation of new grains. Indeed, 
the actual nucleation rate can be defined as 

(8) 

Concerning the limitations of the KJMA theory, it also holds in the case of anisotropic 
growth provided that the grains have a convex shape and are aligned in parallel.— Moreover, 
the KJMA theory provides a good approximation when the anisotropy is moderate or for 
soft impingement.— However, KJMA theory fails when nucleation is non-random,— when 
growth is anisotropic,— i^*^ when growth stops before crystallization is complete^ and when 
the incubation time is not negligible.— 

III. STATISTICAL PROPERTIES OF THE KJMA CELL SIZE DISTRIBUTION 

The cell size distribution is characterized by its PDF, f{s), the probability that a cell has 
a size between s and s + ds. From its definition it is obvious that f{s) must be normalized 



oo 



f{s)ds = 1. (9) 







To analyze the properties of the cell size distribution, we will consider the contribution of 
the crystals formed at a time r over a time interval dr (r-crystals). We will call the cell 
size distribution of the r-crystals the r-distribution. Accordingly, we define the PDF of the 
r-crystals, /r(s), as the probability that a r-crystal has a volume between s and s + ds. 
From the definition of fris), it is also apparent that /t(s) must be normalized: 



POO 

/ fAs)ds = l. (10) 
Jo 



f{s) is simply the addition of the contributions of the r-crystals over the time interval in 
which their nucleation takes place 

m = ^° .o^vv ,y • (11) 

Jo Ia{T)dT 

Note that the denominator in Eq. ( fTTl) ensures that f{s) is normalized if all /r(s) are 
normalized. 

In the following sections, we will present the expected grain size and the variance of the cell 
size distribution and their relationship with the equivalent parameters of the r-distributions. 

A. Expected grain size 

It is well known that the expected grain size, E, is the inverse of the final grain density: 

E = (^j^ h{T)dr^ . (12) 

Likewise, the expected value, Er of a r— distribution is simply the final average grain size 
of a r-crystal normalized to the total volume: 

/oo 
d,v{z,T)dz. (13) 

Introducing Eq. (jl]) into Eq. ( |T3i) leads to 

1 1"^ 

where the extended average growth rate is given by 

d-,VeAz,T) = Dgnr{z,Tf-^G{z). (15) 

Note than once the evolution of the transformed fraction, X(t), is known - i.e. the 
solution of Eq. ([3]) - the calculation of Er is straightforward. 



Besides, the final space fraction occupied by tlie r-crystals, X^, can be calculated from 
the integration over the entire space of the probability that a point P in the space belongs to 
a tau crystal nucleated at O. Since the system is homogeneous and isotropic, this probability 
only depends on the distance h between O and P. Therefore, 

Xr = Iir)jp,{0,P)dVp = 

= DgoI{r) / Pr{b)b''-'db, (16) 

where Prib) is the probability that a point P, separated by a distance b from the nucleus 
O, belongs to the crystal nucleated at O. To calculate Prib), we will use the same approach 
that Kolmogorov used for the deduction of Eq. ([3]). Since nucleation is Poissonian, Pr{b) is 
given by the probability that no nucleus is formed that could transform P before O does so. 
P would be transformed by O at the moment tb- 

b = r{tb,T) = j^' G{z)dz. (17) 

Thus, the nuclei formed at z that could transform P before O does so, are located in a 
/^-sphere of radius r{tb, z) around P. Therefore, according to Eq. (j3]), -Pr(^) is given by 



Prib) = exp 



go I I{z)r{tb,z)^dz 







(18) 



The previous integral spans the time interval [0, tb] since no nucleus formed after tb could 
transform P. Comparison of Eq. (JTj) with Eq. flTSl) gives 

Pr{b) = exp [~X,,{tb)] = 1 - X{tb). (19) 

Finally, if we introduce the value of Prib) given by Eq. (ITSl) into Eq. ( |T6i) and we change 
the variable b by tb, we obtain 

/oo 
[l-X{tb)]r{tb,T)''-'Gitb)dtb. (20) 

Alternatively, the expected value, Er is the ratio between the space fraction occupied by 
the r-crystals and the density of r-crystals: 
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As expected, substitution of Eqs. (120|) and ([8]) into Eq. (12T|) delivers Eq. (fT^ . Moreover, 
the integration of Xt- over the whole time interval where nucleation takes place gives the 
total transformed fraction, 1: 



Ia{r)ErdT = / XrdT = 1. (22) 
Jo 

We will end this subsection verifying that the value of E evaluated from the r-distributions 
coincides with the value given at the beginning of this subsection [Eq. ( fT2i) ]: 

E^ / sfis)ds = ^° ;y ^ ^ = 

Jo Ia{T)dT 

(23) 





j^h{T)ErdT_ 1 



h{r)dT Ia{T)dT 
B. Variance of the grain size distribution 

To determine the variance we will adapt the development of Gilbert^ for a Poisson-Voronoi 
tessellation to our case. First, we define a new PDF, f*{s), as the PDF of the crystals that 
contain a given arbitrary point O: i.e., if we pick an arbitrary point O, f*{s) is the probability 
that a crystal has a size between s and s + ds and contains the point O. Accordingly, f*{s) 
is proportional to f{s) and to s, because a large crystal has a proportionally greater chance 
of containing the point O. Therefore, 

ris) = (24) 
The constant of proportionality, -E"^, has been deduced by imposing normalization: 

ris)ds = l. (25) 
From the definition of f*{s) it can be easily proved that 

POO 

var= / {s- EYf{s)ds = E*E ~ E^ (26) 



Jo 

where E* is the expected value of f*{s). 

Once E* is known, the calculation of the variance is simple. To obtain E* we first analyze 
the contribution of the r-crystals. To do so, we define f*{s) as the PDF of the r-crystals 
that contain a given arbitrary point O: i.e., if we pick an arbitrary point O, f*{s) is the 
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probability that a r-crystal has a volume between s and s + ds and contains the point O. 
Accordingly, f*{s) is proportional to /t(s) and to s: 



f:{s)^Uis)s. (27) 

On the other hand, the integration of f*{s) over all possible volumes is the probability 
that an arbitrary point O belongs to a r-crystal. This probability is the fraction of the space 
occupied by the r-crystals Xr'. 

POO 

/ f:is)ds = X^, (28) 

taking into account that, 

poo 

Er= / Sfr{s)ds, (29) 

and combining Eqs. ([27]), ([28]) and ([2ll) we obtain f*{s): 



f:{s) = IaiT)sMs). (30) 
Then the expected value of /*, E*, is 

It can be easily verified that E* is related with E* through 

E* = = r X,E*dT. (32) 



f*{s)ds Jo 

Therefore the contribution of the r-crystal to the expected value E* is X^-E*. In addition, 
this contribution is the integration over the entire space of the probability that a differential 
volume around a point P in the space belongs to the same r-crystal as O. Since the system 
is homogeneous and isotropic, this probability only depends on the distance h between O 
and P, 

/■oo 

X,E; = Dgn P;(6)6^"M6, (33) 
Jo 

where P*{b) is the probability that two points, O and P, separated by a distance b belong 
to the same r-crystal, 

P;(6) = /(r) J P:{b,Q)dVQ. (34) 

The integration domain covers the entire space, dVg is the D-volume differential around 
a point Q, I{T)dVQdT is the probability that a nucleus is formed at Q at the time r during 
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FIG. 1: Schematic drawing of the calculation of P*{b,Q). 



the time interval dr and P*{b, Q) is the probabihty that both O and P belong to the same 
crystal nucleated at Q (see Fig. [1]). For D = 2, dVq = 2rodrod0o while for D = 3, 
dVq = 2nrQsm{9o)drod6o in polar coordinates. 

P*{b, Q) is given by the probability that no nucleus is formed that could transform O or 
P before Q does, then 



P:{b, Q) =exp I - [ ^ ° I{z)r{to, z)''dz + 

+ [ I{z)r{tp,z)^dz- [ I{z)—dz]}, (35a) 
Jo Jo 9d -iJ 

rl = rl + b^ - 2rob cos Oq, (35b) 

i^x = ^(tx,^) = / G{z)dz, ioYx = 0,P (35c) 

ro + rp-b^ ^ r' ^^^^^^^ ^35^^ 



2 

where r{to, z) and r{tp, z) are the minimum distance between O, P and a nucleus created at 
the time z, so that the nucleus would not transform O and P respectively (see Fig. [1]). Vj is 
the volume intersection between two D-spheres of radius r{to, z) and r{tp, z) centered at O 
and P, respectively (gray region in Fig. [T]). The subtraction of the term Vj is in accordance 
with the fact that it has been accounted twice in the first and second integrals in Eq. fl35al) . 
For a particular set of values of the integration variables ro and 6o, fp is evaluated from 
Eq. (135bp while to and tp are defined by Eq. fl35cl) and t' is defined by Eq. ( ]35dl) . Note 
that O and P are transformed by Q at the times to and tp, respectively. Thus, any nucleus 
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formed after to and tp could not transform O or P, respectively, so the two first integrals 
in Eq. fl35al) span the time interval [0,to] and [0,tp], respectively. Additionally, it can be 
easily verified that if 2; > t', then the intersection between the D-spheres is null. Therefore, 
the last integral in Eq. fl35ap spans the time interval [0,t']. 

Finally, Eq. is simplified by substitution of Eq. ([7]) in the first and second integrals 
in Eq. fl35aD 

exp [-Xexito) - Xexitp)] 



p:ib,Q) 



exp — I{z)Vidz 
[l-X{to)] [1-X{tp 



exp 



-j;;iiz)Vrdz 



(36) 



It can be easily proved that the variance of the r-distributions, var,-, is given by 



var, = E^Er - {Erf. 



(37) 



Finally, we will check if the variance of the distribution determined from the decomposi- 
tion of f{s) into r-PDF, Eq. ( ITTl) . gives the expected result, Eq. (1261) : 



var 



Ia{T)dT 



E 
E 



/oo / poo \ 
n sf:is)dsjdT 



-E^ = 
-E^ 



Erdr 



-E^ = E*E - E\ 



(38) 



At this point, we would like to point out that the results obtained so far are exact and 
general, i.e., we have not made any assumption concerning f{s) and /r(s). 



IV. APPROXIMATE VARIANCE 



According to our previous analysis, the exact calculation of the variance is reduced to the 
calculation of the parameters E and E* in Eq. fl26l) . While the calculation of E is straight- 
forward, the evaluation E* is more cumbersome. Indeed, when compared to Monte-Carlo 
algorithms,— its numerical calculation is more complex without representing any significant 
reduction in computing time. That is because there are several integrals nested and, in 
particular, the calculation of the intersection volume Vj is complex. When r{to,z) ^ 6 or 
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b ^ r(to,z), Vj tends towards being a D-sphere of radius r(to,z) and 0, respectively. On 
the other hand, when r{to,z) ~ r{tp,z) the shape of Vj roughly approaches a D-sphere. 
Since the width of Vj (see Fig. [T]) is r(to, z) +r(tp,z) —b, we approximate V/ by a D-sphere 
of diameter r{to, z) + r{tp, z) — b: 

y^^^(rJfo,^)+Z^^^\\ (39) 



2 

It is worth noting that the previous approximation also works for the limiting cases 
r{to,z) ^ b and b ^ r{to,z). Furthermore, for the one-dimensional case it can be easily 
verified that Eq. ( l39i) is exact (in Appendix A we derive P* for D = 1). Finally, the 
approximate solution (from here on approximation I) is obtained by substitution of Eqs. 
and ([39D into Eq. 

exp [-Xexito) - Xex{tp)] 



exp[-X,,(t')] 
[l-Xito)][l-Xitp)] ^^^^ 



[l-Xif) 

Therefore, the calculation of P*{b,Q) is simple provided that the evolution of the trans- 
formed fraction, X{t), is known. Analytical exact solutions for X{t) are restricted to three 
particular situations under isothermal conditions: time- independent growth and nucleation 
rates, time-independent growth rate and nucleation rate proportional to a power of time,— 
and site saturated nucleation. A quasi-exact solution of the KJMA model has recently been 
obtained under continuous heating conditions.— Moreover, there are numerical methods 
which allow a simple and fast calculation of X{t) for an arbitrary time dependence of the 
nucleation and growth rates.— 

We have analyzed the distribution emerging from solid phase crystallization under 
isochronal heating conditions, i.e. heating at a constant rate, to check the accuracy of ap- 
proximation I, Eq. fHOj) . in the case of time dependent nucleation and growth rates. To work 
with realistic parameters we have taken those of amorphous silicon crystallization,— 1^>^ in 
which the nucleation and growth rates are described by an Arrhenius temperature depen- 
dence 

/ = /oexp (^-^^^ and 6" = Go exp (^-^^^ . (41) 

where T is the temperature in Kelvin and kp is the Boltzmann constant. In Table [T] we 
summarize the corresponding parameters. When the temperature is raised at a constant 
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TABLE I: Experimental nucleation and growth rates of amorphous sihcon (Ref. 
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3S 



Nucleation 



Growth 



Activation energy, E^- 
Preexponential term, Iq 
Activation energy, Eq 
Preexponential term. Go 



5.3 eV 

1.7 X lO'^^s-im"^ 
3.1 eV 

2.1 X lO^s^^m 



rate /?, the nucleation and growth rates become time dependent through Eq. ( 14T|) . Under 
those conditions, the kinetics is correctly described by the KJMA theorj^i^ and there is 
good agreement between experiment and theoretical predictions.—*^ For the calculation of 
ution of the transformed fraction, we have used the quasi-exact solution described 



the evo 

in Ref. |55|. The numerical evaluation of the integrals has been performed by means of an 
extended midpoint algorithm.— To confirm that the observed discrepancies are not related 
to numerical inaccuracies, we have performed several calculations with consecutive smaller 
integration steps. Moreover, for the numerical integration over a semi-infinite interval, we 
have imposed a minimum relative error of 10~^. To check the accuracy of the numerical 
calculation, we have calculated the integral of Xr over the interval [0, oo) and have compared 
them to its predicted value, Eq. fl22|) . Calculations that exhibit discrepancies larger than 
10~^ were rejected. 

As is apparent from Fig. [H the approximation of Vj by a D-sphere of diameter equal to 
its width, Eq. fl39l) . results in an underestimation of Vj, which leads to an undervaluation 
of E* and of var^. The latter conclusion can be verified in Fig. [21 where the evolution of 
E* and var,- with r is shown. Although approximation I gives an accurate value of E*, the 
approximate value of var,- shows a significant deviation from the exact value. The reason is 
that in the evaluation of var^, Eq. flTTI) . both terms in the difference have similar values. The 
same happens to the values of the variance and E*; the exact and approximate values of E* 
are 3.93 and 3.69 respectively, while the exact and approximate variances are 3.56 and 3.22, 
respectively. (Space has been normalized to the space scaling factor— {G{Tp) / I{Tp)Y^'^, 
where Tp is the peak temperature.) Despite the significant discrepancy between the exact 
and approximate values of var^, they have a nearly parallel evolution with r. This result 
is general and is related to the very similar dependency of the approximate and exact Vj 
on the integration parameters. Therefore, the accuracy of approximation I can be analyzed 
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FIG. 2: E* and var,- for three-dimensional growth and isochronal heating of 40 K/min. Time and 
space have been normalized according to the time and space scaling factors, {I{Tp)G{Tpf)-^/^ 
and {G{Tp)/I{Tp)y/'^, where Tp is the peak temperature, i.e., the temperature at which the 



transformation rate is maximum (see Ref. 
(red dashed line) values are compared. 



35l ). The exact (black solid line) and the approximate 
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En/Eg 



FIG. 3: Exact (black triangles) and approximate (red circles) values of varo and their ratio (solid 
blue squares) as a function of the ratio between nucleation and growth activation energies. 



through the relation between the exact and approximate values of vaiQ. To cover a wide 
range of distributions, we will recall the results given in Ref. Issj. In this work it was shown 
that the shape of the grain size PDF was practically insensitive to the heating rate, but 
it depends mainly on the ratio En/Eq, i.e., the relative evolution of the nucleation and 
growth rates with time. The limit E^/Eg — > corresponds to site saturated nucleation 
while En/ Eg = 1 coincides with the isothermal case. 
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In Fig. [3] we have plotted the exact and approximate values of varo as well as their ratio. 
(At this point it is worth recalling that, according to Eq. fHTl) . a relation of one order of 
magnitude between the activation energies E^q and Eq would result in a huge difference in 
the relative time evolution between the nucleation and growth rates.) First, we can easily 
verify that varo (and in general the variance) decreases with Ej^i/Eg- This means that the 
distributions become broader as En /Eg increases. Indeed, when Eg ^ -Eat, during the 
first stages of the transformation, nucleation dominates over growth. Most of the nuclei are 
formed at the beginning and they grow at a slow rate. Thus, the average grain size and its 
variance diminishes when En /Eg diminishes. In contrast, when En ^ Eg, during the first 
stages of crystallization, growth dominates and the nucleation rate increases progressively as 
crystallization proceeds. Since the time left for growth is less for the nuclei that appear later, 
the density of small grains will be higher than that of larger grains. So the average grain size 
and the distribution variance increase with En /Eg- On the other hand, despite the large 
variation of varo, the total variation of the ratio between the exact and the approximate 
varo is very smooth - from 1.96 to 2.12. A similar behavior has been observed for the 2D- 
case where this rate evolves from 1.28 at En /Eg = 10 to 1.35 at En /Eg = 0. Hence, the 
deviations of the approximate value of varo from the exact value remain practically constant. 
This result is due to approximation /, which is based on a geometrical approach that is fairly 
insensitive to the relation between nucleation and growth rates. 

Since the ratio between the exact and the approximate Yax-r is nearly constant, we can 
obtain a significantly more accurate approximate value for var^ by simply multiplying it 
by the corresponding proportionality constant. This constant only depends on the growth 
dimensionality. We have chosen the values of 2.07 and 1.32 for the 3D and 2D cases, 
respectively. These values correspond to En/ Eg = 1, i.e., they correspond to the isothermal 
case with Eg and En constant in time. With this correction (from now on approximation 
II), the relative error in the calculation of the variance diminishes to less than 2% and 
for En /Eg > 1 the relative error is less than 0.2%. In Fig. H] we have plotted the exact 
and the approximate values of the distribution variance and E* with respect to the ratio 
En/ Eg- The exact and the approximate value obtained from approximation II of the 
variance and E* exhibit excellent agreement; the values overlap in such a way that they are 
nearly indistinguishable. Concerning the values obtained from approximation I, it is worth 
noting that despite the significant error related to the calculation of varo [Fig- E] and of var,- 
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En/Eg 

FIG. 4: Exact (black triangles), approximate (red circles) and corrected approximate (blue squares) 
values of the variance and E* and as a function of the ratio between the nucleation and growth 
activation energies. 

in general, the inaccuracies in the evaluation of the variance and E* are significantly smaller. 
The reason is that both parameters depend exclusively on and E [Eqs. (1261) and fl32|) ]. 
From Fig. [2] it is clear that the error related to E^- is significantly smaller than the error in 
the evaluation of var^, while the calculation of E is exact. 

From now on, when will always use approximation II in the calculation of the approximate 
values of E* and E*. 

V. APPROXIMATE CELL SIZE PROBABILITY DENSITY FUNCTION 

One application of the preceding analysis of the statistical properties of grain size dis- 
tribution is the derivation of a PDF. If we choose a set of /r(s) such that their expected 
value coincides with the result of Eq. f|T4|) and their variance is equal to the value given by 
Eq. (!37l) . then the variance and expected values of the PDF obtained from Eq. ( fTTl) will be 
exact, i.e., the PDF obtained from Eq. ( ITTl) will have the same variance and expected values 
as the actual PDF. Indeed, Pineda et al.,— apply this approach in the case of tessellations 
generated by random nucleation processes where the growth rate was assumed to be con- 
stant. The agreement between their approximate PDF and Monte-Carlo simulations was 
remarkable. However, they did not notice that the expected value and the variance of their 
approximate PDF were exact. Concerning the particular choice of the /r(s) functions, we 
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FIG. 5: Grain size distribution for three-dimensional growth and isochronal heating. Comparison 
between the PDFs obtained from the exact calculation of E* (black solid line), from the approxi- 
mate calculation of E* (red dashed line) and from Monte-Carlo simulation (empty squares). 

will consider two cases: gamma and Gaussian distributions. 
A. Gamma distribution 

Given that in a r-distribution the nucleation events are simultaneous and the r-nuclei are 
randomly distributed, we can assume that /t(s) is similar to the PDF resulting from a process 
of site saturated nucleation, i.e. that the PDF is that of a Poisson-Voronoi tessellation. 
As explained in the introduction, the gamma distribution is the exact PDF for the one- 
dimensional case while it provides a very accurate result for the two- and three-dimensional 
cases: 



Since E^- is already the expected value, the problem comes down to the determination of 
the exponent z/,-. Indeed, can be calculated from the following property of the gamma 
distributions 



To check the accuracy of the PDF we compare them to some Monte-Carlo simulations. 
The Monte-Carlo algorithm^''' consists in dividing the space into a cubic lattice. Cells are 
assigned to nuclei randomly. The nucleation time of each nucleus is precisely calculated from 





(42) 



var.^ 
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FIG. 6: Grain size distribution for two-dimensional growth and isochronal heating. Comparison be- 
tween the PDFs obtained from the exact calculation of E* (black solid line) , from the approximate 
calculation of E* (red dashed line) and from Monte-Carlo simulation (empty squares). 

the nucleation rate and it is recorded to evaluate the exact evolution of the grain growth. 
Finally, each cell is assigned to the nucleus that first reaches this cell. The evolution of 
the grain growth transformation is checked whenever the grain growth is equal to the size 
of a cell in order to avoid incorrect cell assignation due to shielding effects. In particular 
we consider the crystallization of amorphous silicon under isochronal heating at 40 K/min 
(Table [T]). The results of the calculations for three- and two-dimensional growth are given in 
Figs. ([5]) and IQ, respectively. We have calculated the PDF using the exact, Eq. (136|) . and 
approximate, Eq. ( I40l) . values of E*. We will refer to these PDFs as the approximate III and 
approximate IV PDFs, respectively. The validity of the selection of a gamma distribution 
for the /t(s) functions is confirmed by the good agreement between the calculated PDF 
and the Monte-Carlo simulations. The excellent agreement between the approximate III 
and the approximate IV PDFs is also noteworthy - the relative difference is less than 0.1% 
for s/E > 0.05. Therefore, the approximate calculation of E* is useful to obtain a simple 
and accurate PDF for a KJMA tessellation. However, the complexity and the computing 
time required for their evaluation is significantly different. For instance, the calculation 
of the approximate III PDF typically takes more than thirty times the time required for 
the calculation of the approximate IV PDF. The reason for this significant simplification 
and reduction in computing time is that the approximate IV PDF saves us from having to 
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FIG. 7: Grain size distribution for three-dimensional growth and site saturated nucleation. Com- 
parison between the PDFs obtained from the exact calculation of E* (black solid line), from the 
approximate calculation of E* (red dashed line) and from Monte-Carlo simulation (empty squares). 

evaluate the inner and more complex integral of Eq. (1361) . 

To confirm this last conclusion, we have calculated the PDF for the case of site saturated 
nucleation. This case corresponds to En /Eg = in Figs. ([3]) and (jlj) and is the case that 
exhibits the greatest discrepancy between the exact and approximate values of the variance 
and E*. Therefore, it should give the worst agreement between approximate III and the 
approximate IV PDFs. As is apparent from Fig. ([7]) here again the agreement between the 
PDFs obtained from the exact calculation of E*, from the approximate calculation of E* and 
from Monte-Carlo simulation is excellent. Only small deviations of the approximate III from 
the approximate IV PDF are distinguishable for s ^ E. Finally, for the one-dimensional 
case and for site saturated nucleation, both the approximate III and approximate IV PDFs 
turn into the exact PDF (see Appendix B). 

B. Gaussian distribution 

When analyzing the crystallization morphology, it is often better to use the grain radius 
distribution instead of the grain size distribution. The grain radius, r, of a grain of size 
s, is defined as the radius that will have a D-sphere of volume s. For instance, for three- 
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dimensional growth 



Given the grain size PDF, the grain radius PDF, g{r), can be easily derived; for three- 
dimensional growth, 

g{r) = Anr'^f i-nr^] . (44) 



^3 

On the other hand, for site saturated nucleation and three-dimensional growth it has 
been shown that g{r) is accurately described by a Gaussian distribution^ 

where Eg is the expected value of the grain radius PDF and a is the standard deviation 



cr 



yvaf^. (46) 



where var^ is the variance of the grain radius distribution. The fitted parameters were 
actually Eg = and a = 0.0883/nQ^^, where no is the nuclei density. 

Conversely, g{r) can be determined by means of Eq. (jH]) where f{s) is a gamma distri- 
bution with u = 5.581.— Under this approach, it can easily be proved that 

1/3 



47r z/ / T(u) 



3 E\ 



2/3 



r(i/ + 2/3) r(z/ + i/3) 



(47) 



where E = I/hq is the expected value of the grain size distribution (see Appendix B). 



1 /3 

By substituting u = 5.581 into the previous relations we obtain Eg = 0.608/r;,Q , var. 
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7.6510-7^0^^ and a = 0.0874/r^y^ which are in good agreement with the parameters 
obtained from the Gaussian fit. 

Therefore, for a KJMA tessellation we can obtain g{r) as the superposition of grir) 
PDFs. In contrast with the previous subsection, we will now assume that grir) are Gaussian 
distributions: 

. ^ /o°°/a(r)^^(r)c;r 
Jo lairjdr 

To evaluate the gri^) PDFs we need to know their expected value, Eg^r, and their vari- 
ance, var^ These parameters can be easily derived from the statistical parameters of the 
r-distributions, E^ and var,-, by means of Eq. (j47j) . In Fig. ([8]) we have plotted the grain 

20 




FIG. 8: Grain radius distribution for three-dimensional growth and isochronal heating. Comparison 
between the PDF obtained from Eq. ()44p (black solid line), Eq. (j48p (red dashed line) and Monte- 
Carlo simulation (empty squares). 

radius distribution obtained, from the grain size distribution [Eq. (jll])] (where f{s) is the 
superposition of gamma distributions), as a direct superposition of grir) Gaussian PDFs 
[Eq. fHH]) ] and from Monte-Carlo simulations. This distribution corresponds to the crystal- 
lization of amorphous silicon under isochronal heating at 40 K/min (Table [T]). The good 
agreement between the PDF calculated from Eqs. (jUj) and fHS]) confirms that the Qrir) are 
correctly described by a Gaussian distribution. Finally, both approaches show good agree- 
ment with Monte-Carlo simulations, i.e., both approaches are useful for describing the grain 
radius PDF. 

VI. CONCLUSIONS 

This paper deals with a subject of interest in many scientific areas, namely the cell- 
size distribution of space tessellations that emerge from first-order phase transformations 
ruled by nucleation and growth of the new stable phase. No restrictions are imposed on 
the time dependency of the nucleation and growth rates, and the validity of our results is 
limited to transformations that obey the premises of the Kolmogorov- Johnson- Mehl-Avrami 
model. We have derived some important statistical properties such as the expected value 
and the variance. The approach used is an extension of the work of Gilbert .i Additionally, 
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we have developed a significantly simpler relation for the calculation of the variance. The 
discrepancies between the exact and approximate variances are less than 2%. 

Like Pineda et al.—- we have derived an approximate grain size PDF as the superposition 
of gamma distributions. We have proved that the expected value and variance derived from 
this approximate grain size PDF are exact. Moreover, we have checked its accuracy against 
Monte-Carlo simulations for a system undergoing a crystallization under isochronal heating 
conditions. The results show a remarkably good agreement between the approximate PDF 
and the Monte-Carlo simulations. Finally, we have shown that the grain radius PDF can be 
expressed as the superposition of Gaussian distributions. 

APPENDIX A: ONE-DIMENSIONAL PDF 

For the one dimensional case, the extended transformed fraction is 



{z)dz 



dr, (Al) 



and the expected value Er becomes 



Er = ^TTT / [1 - ^i^)] G{z)dz. (A2) 



l-X(r) 

With regard to the calculation of P*{b), we have split the entire space into three regions. 
The first corresponds to Q located at the left side of O, in this case rp = ro + b and 

P:{b,Q)=exp[-X,,{tp)]. (A3) 

The second region corresponds to Q located between O and P, then rp + vq = b and 

P;(6, Q) = exp[-Xe.(to) - X,^{tp) + Xe.(r)]. (A4) 

And the third region corresponds to Q located at the right side of P, tq = rp + b and 

P;(6,g)=exp[-Xe,(to)]. (A5) 

Finally, combining Eqs. (\A3\i . flA4|) and flASP with (|34|) we obtain 

P;(6) = Iir){2j^^[l-X{to)]dro + 
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ro 



to 



G{z)dz, 
tp 



Once P*{b) is known, Eq. 



h — To = J G{z)dz. 
combined with Eq. fl571) delivers var^ 



(A6) 



APPENDIX B: SITE SATURATED NUCLEATION 

When nucleation is completed prior to crystal growth, the nucleation rate can be approx- 
imated to I{t) ^ no6{t) where no is the density of nuclei and 6 is the Dirac delta function. 
In this case, the extended transformed fraction becomes 

D 



Xexit) = nogo 



G{z)dz 



Then 



4(r) = no[l-X(r)]5(r) 

and the PDF, Eq. ( fTTj) . is reduced to 

_ j^n,[l-X{rmr)U{s)dr _ 



(Bl) 



(B2) 



(B3) 



Therefore, in this case we only need to calculate Eq and varo. Concerning Eq, Eq. f|T^ . it 
can be easily proved that it is simply 

1 



Eq 



no Jo 



exp[~-Xex]dXex = l/riQ. 



(B4) 



Indeed, according to Eq. fIBSp Eq = E, and the expected value of E is I/uq [Eq. f|T2|) ]. 
Moreover, the fraction of space occupied by a r-crystal, X^-, is reduced to 5(r), as expected. 

With respect to varg, its value is determined by E^ and Eq by means of Eq. flTTI) . Eq 
is given by Eq. (1331) . We therefore need to evaluate PQ{b), the value of which depends on 
which relation for PQ{b, Q) we use: the exact one Eq. (136!) or the approximate one Eq. (HOl) . 
We will evaluate Poib) for the one-dimensional case because in this case the approximate 
solution coincides with the exact one. Specifically, substitution of Eq. ( IBll) into Eq. (1A6I1 
leads to 



P*ib) = no6{T) 



POO 

'I 

2nob , 



O 



no5(r)e-'"«'' {I/hq + b) . 



(B5) 
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Then, from Eqs. (133!) and (IBSP 

° 2no 

And finally, from Eqs. (1371) and ( lB6p . 

varo = (2rao)~ 

If we choose a gamma distribution for the calculation of the cell-size PDF, from Eq. ( l42l) we 
obtain uq = 2. Thus for site-saturated nucleation and one-dimensional growth we obtain a 
gamma distribution with u = 2 and E = I/uq which agrees with the exact solution.— 
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